Screening effects in flow through rough channels 



J. S. Andrade Jr. 1 ' 2 ' 3 , A. D. Araujo 1 , M. Filoche 2 ' 3 , and B. Sapoval 2 ' 3 
1 Departamento de Fisica, Universidade Federal do Ceard, 
60451-970 Fortaleza, Ceard, Brazil 
2 Centre de Mathematiques et de leurs Applications, 
Ecole Normale Superieure de Cachan, 
94235 Cachan, France 
3 Laboratoire de Physique de la Matiere Condensee, 
Ecole Polytechnique, 
91128 Palais eau, France 
(Dated: February 4, 2008) 

A surprising similarity is found between the distribution of hydrodynamic stress on the wall 
of an irregular channel and the distribution of flux from a purely Laplacian field on the same 
geometry. This finding is a direct outcome from numerical simulations of the Navier- Stokes equations 
for flow at low Reynolds numbers in two-dimensional channels with rough walls presenting either 
deterministic or random self-similar geometries. For high Reynolds numbers, when inertial effects 
become relevant, the distribution of wall stresses on deterministic and random fractal rough channels 
becomes substantially dependent on the microscopic details of the walls geometry. In addition, we 
find that, while the permeability of the random channel follows the usual decrease with Reynolds, 
our results indicate an unexpected permeability increase for the deterministic case, i.e., "the rougher 
the better". We show that this complex behavior is closely related with the presence and relative 
intensity of recirculation zones in the reentrant regions of the rough channel. 

PACS numbers: PACS numbers: 05.45.Df, 41.20.Cv, 47.53.+n 



Partial differential equations are basic in the mathe- 
matical formulation of physical problems. The Laplace 
equation, for example, is known for its relevance in many 
fields, namely electrostatics, heat transport, heteroge- 
neous catalysis, and electrochemistry. Although rather 
simple in form, the Laplace equation can have highly 
non-trivial solutions, specially if the boundary of the sys- 
tem represents the surface of an irregular object. The 
role of this particular complexity has been a theme of 
constant research interest with recent important devel- 
opments [l|, |3, Isj. Much more complex, however, are 
the Navier- Stokes equations for the description of hydro- 
dynamic flow in irregular geometries. The aim of the 
present work is first to reveal a surprising analogy be- 
tween the properties of the solutions of 2d Laplace and 
Navier-Stokes equations when flow at low Reynolds and 
proper boundary conditions are imposed on the same ge- 
ometry. In a second step, we present results obtained at 
higher Reynolds numbers and for distinct types of surface 
geometry. 

The situations we compare and find to be quantita- 
tively similar are displayed in Fig. 1. Figure la pictures 
the simplest Laplacian problem, namely that of a capac- 
itor with a pre- fractal electrode. The complex feature 
here is the distribution of charge on the irregular elec- 
trode or, in mathematical terms, the distribution of the 
harmonic measure. The recent research on this field has 
been mainly dedicated to the application of Laplacian 
transport towards and across irregular interfaces p], 0, j5| . 
The system depicted in Fig. lb corresponds to Stokes flow 
in a symmetric rough channel with the same geometry as 
in Fig. la. The hydrodynamic quantity which is found 
to be distributed similarly to the harmonic measure in 
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FIG. 1: The two different problems with similar solutions. In 
(a) we show a capacitor with an irregular electrode. The local 
charge is obtained from the numerical solution of the Laplace 
equation with Dirichlet boundary conditions. In (b) we see 
the analogous problem of flow at low Reynolds number. The 
stress parallel to the wall channel can be calculated from the 
solution of the continuity and Stokes (Re — 0) equations. 



Fig. la is the viscous stress along the channel boundary. 

To obtain the potential field for the problem in Fig. la, 
one must compute the solution of Laplace equation with a 
potential V = 1 on the counter electrode and zero poten- 
tial on the irregular electrode (Dirichlet boundary con- 
dition). From that, it is then possible to calculate the 
charge a l L on each elementary unit i of the wall and its 
corresponding normalized counterpart, (j) l L = o- l L /^(j 3 L . 
As previously shown 6], the distribution of 4>l along 
the irregular electrode is strongly nonuniform as a conse- 
quence of screening effects, being in fact characterized as 
multifractal. This is a typical situation where the deep 
regions of the irregular surface support only a very small 
fraction of the total charge, as opposed to the more ex- 
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posed parts. 

The hydrodynamic question posed here has been trig- 
gered by the idea that the design of flowing systems 
should include the influence of the surface geometry as 
a possibility for optimal performance. For this, we in- 
vestigate the flow in a duct of length L and width h 
whose delimiting walls are identical pre-fractal interfaces 
with the geometry of a square Koch curve (SKC) Q- 
The mathematical description for the fluid mechanics in 
this channel is based on the assumptions that we have 
a continuum, Newtonian and incompressible fluid flow- 
ing under steady state conditions. The relevant physical 
properties of the fluid are the density p and the viscos- 
ity jjL. In our simulations, we consider non-slip bound- 
ary conditions at the entire solid- fluid interface. In ad- 
dition, gradientless boundary conditions are assumed at 
the exit x = L, whereas a parabolic velocity profile is 
imposed at the inlet of the channel. The numerical so- 
lution for the velocity and pressure fields in the rough 
channel is obtained through discretization by means of 
the control volume finite- difference technique 7] . For our 
complex geometry, this problem is solved using a struc- 
tured mesh based on quadrangular grid elements. For 
example, in the case of the channel with walls that are 
fourth-generation SKC curves, a mesh of approximately 
a million elements adapted to the geometry of the inter- 
face generates satisfactory results when compared with 
numerical meshes of higher resolution. 

The Reynolds number is defined here as Re = pVhj '/i, 
where V is the average velocity at the inlet. In Fig. 2a 
we show the velocity vector field at low Re located at the 
self-similar reentrant zones that constitutes the rough- 
ness of the irregular channel. Indeed, as depicted in 
Fig. 2b, by rescaling the magnitude of the velocity vectors 
at the details of the roughness wall, we can observe fluid 
layers in the form of consecutive eddies. Although much 
less intense than the mainstream flow, these recirculating 
structures are located deeper in the system, and therefore 
experience closer the landscape of the solid-fluid inter- 
face. More precisely, viscous momentum is transmitted 
laterally from the mainstream flow and across successive 
laminae of fluid to induce vortices inside the fractal cav- 
ity. These vortices will then generate other vortices of 
smaller sizes whose intensities fall off in geometric pro- 
gression jslEJ. 

Once the velocity and pressure fields are obtained for 
the flow in the rough channel, we can compute the nor- 
malized stress <j) l s at each elementary unit i of the wall, 

(j) % s = Tg/J2 T si wnere the sum is over the total num- 
ber L p of perimeter elements, the magnitude of the local 
stress is given by rs = \dv\\/dn\, the derivative is calcu- 
lated at the wall element, v\\ is the local component of 
the velocity that is parallel to the wall element, and n is 
the local normal coordinate. The semi-log plot in Fig. 3a 
shows that the spatial distribution of normalized stresses 
at the interface is highly heterogeneous, with numerical 
values in a range that covers more than five orders of 
magnitude. Also shown in Fig. 3a is the variation along 




FIG. 2: (a) Vortices in the reentrant zones of the upper half 
of the (symmetric) pre-fractal roughness channel. The fluid 
flows steadily from left to right at a low Reynolds number, 
Re = 0.01. (b) Sequence of the smaller eddies at the details 
of the roughness wall shown in (a). These structures can only 
be viewed by rescaling the magnitude of the velocity vectors 
because their intensities fall off in geometric progression. 



the interface of the normalized Laplacian fluxes (j)L cross- 
ing the wall elements of a Laplacian cell with Dirichlet 
boundary conditions. The astonishing similarity between 
these two distributions clearly suggests that the screen- 
ing effect in flow could be reminiscent of the behavior 
of purely Laplacian systems. As shown in Fig. 3b, this 
analogy is numerically confirmed through the very strong 
correlation between local stresses and Laplacian fluxes. 
These measures follow an approximately linear relation- 
ship, namely (ps <\>l- 

A further analogy can be drawn from the discussion of 
the notion of active zone, as the zone which supports 
the majority of the charge or current (for the Lapla- 
cian field) as well as the shear stress here lj. For two- 
dimensional systems subjected to Dirichlet's boundary 
condition, a fundamental step towards the understand- 
ing of the purely Laplacian problem is the mathematical 
theorem given by Makarov 2] . This theorem states that 
the information dimension of the harmonic measure on a 
singly connected interface of arbitrary geometry in D = 2 
is exactly equal to 1. In practical terms, it essentially says 
that, whatever the shape (perimeter) of an interface, the 
size of the region where most of the activity takes place 
is of the order of the overall size L of the system. Here 
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FIG. 3: (a) The (dark) black line is the curvilinear distribu- 
tion of the logarithm of the normalized shear stresses <j)s on 
the interface along one of the two (symmetric) square Koch 
curves corresponding to the channel walls. The (light) red 
line gives the distribution of the logarithm of the normalized 
Laplacian charges <j)L for the analogous electrostatic prob- 
lem. For better visualization, the distribution of (j)L has been 
shifted downwards, (b) Double-logarithmic plot of (j)s versus 
<pL, with the red line indicating their linear relationship. 



we define an active length L a as, 

L a = l/J2(^>s) 2 (l<L a <L p ). (1) 

1 = 1 

If L a is equal to the wall perimeter L p: the entire wall 
works uniformly. The theorem of Makarov indicates that, 
on the opposite, for a purely Laplacian field, L a w L. 
The active length therefore provides an useful index to 
quantify the interplay between the flow and the complex 
geometry of the interface at the local scale. The results 
in Fig. 4a (open circles) show that the value of L a for 
the square Koch curve remains approximately constant 
at L a /L = 0.55 for low and moderate Reynolds numbers. 
This value is consistent with our screening analogy be- 
cause it indicates that the hydrodynamic stress is mainly 
concentrated in a subset of the wall whose size is of the or- 
der of the system size L. Only at higher Re values, when 
inertial forces become comparable to viscous forces, one 
can observe a small increase in L a . The stress becomes 
slightly less localized due to the higher relative inten- 
sities of the vortices inside the deeper reentrant zones, 
when compared with the intensities of the correspondent 
flow structures at low Reynolds conditions. 

The usual approach to describe single-phase fluid flow 
in irregular media (e.g., porous materials and fractures) 
is to characterize the system in terms of a macroscopic 
index, namely the permeability K, which relates the av- 
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FIG. 4: (a) Semi-log plot showing the dependence of the ac- 
tive length L a on the Reynolds number Re for the square 
Koch channel (empty circles) and random Koch channel 
(stars). In (b) we show the semi- log plot of the variation 
of the normalized permeability with Re for the random Koch 
curve (stars) and the second (triangles), third (diamonds) and 
fourth (circles) generations of the square Koch channel. One 
observes that the higher is the generation of the pre-fractal 
SKC, the larger is its permeability. 



erage fluid velocity V with the pressure drop AP mea- 
sured across the system, V = —KAP/jiL. Figure 4b 
(open circles) shows that the permeability of the rough 
channel for low and moderate Re remains essentially con- 
stant at a value that is slightly above but very close to 
the reference permeability of a two-dimensional smooth 
channel, namely (K/K ~ 1), with K = h 2 /12 0. 
Above a transition point at about Re ~ 10, the change 
in permeability reflects the onset of convective effects in 
the flow, and therefore the sensitivity of the system to 
these inertial nonlinearities. Surprisingly, one observes 
that, instead of decreasing with Re (i.e., a behavior that 
is typical of disordered porous media), the permeability 
of the SKC substantially increases. Moreover, as shown 
in Fig. 4b, the higher the generation of the SKC, the 
higher is the permeability of the channel for a fixed value 
of Re above the transition. In other words, "the rougher 
the better" . These results show that the screening effect 
of the hierarchical SKC geometry on the flow can be un- 
derstood in terms of a reduction in the effective non-slip 
solid- fluid interface. In other words, we can imagine that 
each vortex present in a given generation of the SKC is in 
fact replacing one or a set of highly dissipative (non-slip) 
wall elements of previous SKC of lower generations. 

Next we study the fluid flow through a rough channel 
whose walls are composed of 10 successive and distinct 
realizations of the random Koch curve (RKC) of third 
generation. As in the previous case with the determinis- 
tic SKC walls, this task is performed here through finite- 
differences Q, but now with an unstructured mesh of 
triangular grid elements based on a Delaunay network. 
Interestingly, the results shown in Fig. 4a (full circles) 
indicate that the active length L a of the wall stress cal- 



4 



LJL 



- 


j, AA 

]{ i J ( 
Y i V 








1, 

jU WiA/l A 

u 

r 


\ 1 [ . 
- 


/ 


2 3 


4 


5 




r i* 

7 5 


9 10 



50 100 150 200 250 

xll 



FIG. 5: The active length L a for each of the 10 wall subsets 
composing the entire geometry of the random Koch curve 
channel shown at the top. The dashed lines indicate the 
boundaries between consecutive subsets, as numbered below. 
The fluid flows from left to right at Re = 0.1. 

ciliated for the entire irregular interface geometry at low 
values of Re, L a /L ~ 0.46, is not substantially differ- 
ent from the SKC case, where L a /L = 0.55. In Fig. 5 
we show that the ratio L a /L calculated individually for 
each of the 10 wall subsets composing the rough channel 
does not vary significantly from one unit to another. This 
is a rather unexpected behavior, specially if we consider 
the complexity of the different geometries involved (see 
Fig. 5). A similar effect has been observed for the purely 
Laplacian problem in a random geometry . These facts 
strongly reinforce our conjecture that the rationale be- 
hind the screening effect in flow through rough channels 
should be based on a general conceptualization that is 
quite similar to the well established theory for Laplacian 
systems 

As depicted in Fig. 4a, by increasing the Reynolds 
number, the departure from Stokes flow due to convec- 
tion at Re ~ 10 results initially in the decrease of L a /L 
(calculated over the entire surface) down to a minimum 
of approximately 0.25 at Re c± 200. This behavior indi- 
cates the presence of long-range flow correlations imposed 
by inertia among successive wall subsets. More specif- 
ically, due to the randomness of these interface units, 
we can observe either "inward" (e.g., the wall subsets 2, 
3, 4, 5, 6 and 8 in Fig. 4) or "outward" protuberances 
(e.g., the wall subsets 1, 7, 9 and 10 in Fig. 4) composing 



the roughness of the channel. Due to the symmetry of 
the system, the inward elements generate bottlenecks for 
flow. At high values of Re, the effect of inertia is to in- 
duce flow separation lines between the mainstream flow 
at the center of the channel and the flow near the wall, 
that can be as large as the largest distance between two 
consecutive bottlenecks. These exceedingly large "stag- 
nation regions" are responsible for the initial decrease in 
the active length. If we increase even more the Reynolds 
number, the relative intensities of the vortices in these 
regions starts to increase. As in the case of the SKC 
channel, L a /L starts to increase due to a more distribu- 
tion of shear stress near the wall. 

Finally, in Fig. 4b we show the variation with Re of the 
permeability for the RKC channel (full circles). In this 
case, we also observe a transition from linear (constant 
K) to nonlinear behavior that is typical of experiments 
with flow through real porous media and fractures [TH . 
Contrary to the results obtained for the SKC channel, 
however, the value of K calculated for low Re values 
is significantly different and smaller than the reference 
value Kq of the corresponding smooth channel (Kj Ko ~ 
0.6). Once more, this is a consequence of the presence 
of several bottlenecks in the channel, which drastically 
reduce the effective space for flow. At high Reynolds, 
this difference is amplified due to inertial effects. 

In summary, we have investigated the effect of deter- 
ministic as well as random roughness of 2D channels on 
local as well as macroscopic flow properties. Our main 
results are threefold. First, at low Reynolds numbers 
there exists a very close analogy between the spatial dis- 
tribution of the local stress on the rough walls and the 
distribution of charge resulting from the solutions of the 
Laplace equation for the electric potential in the same 
geometry. Second, for a fractal deterministic roughness, 
a surprising increase of the permeability with Reynolds 
is observed. Moreover, this effect is augmented by in- 
creasing the fractal generation of the channel wall so 
that, paradoxically, "the rougher the better". Such re- 
sults could find practical applications in microfluidics and 
help, for example, to understand the enhanced hydrody- 
namical features underlying shark-skin effects. 
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